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SCALING  LAWS  FOR  THE  SPECTRUM  OF  INTERCHANGE 
INSTABILITIES  IN  THE  HIGH  LATITUDE  IONOSPHERE 


1.  INTRODUCTION 

In  the  last  several  years  a  considerable  amount  of  both  experimental 
[Basu  et  al.,  1984;  Basu  et  al,  1988;  Gurnett  et  al.,  1984;  Kintner  et  al., 
1987;  Weimer  et  al.,  1985;  Vickrey  et  al.,  1980;  Vickrey  et  al.,  1986; 
Cerisier  et  al.,  1985;  Baker  et  al.,  1986;  Curtis  et  al.,  1982]  and 
theoretical  [Mitchell  et  al.,  1985;  Lysak  and  Carlson  1981;  Lotko  et  al., 
1987;  Chaturvedi  and  Huba,  1987;  Keskinen  et  al.,  1988]  research  has  been 
directed  to  the  origin,  modeling,  and  general  interpretation  of  plasma 
turbulence  and  structure  in  the  high  latitude  ionosphere  and  magnetosphere. 
It  is  now  known  [for  recent  reviews,  see  Kintner  and  Seyler,  1985;  Temerin 
and  Kintner,  1988;  and  Tsunoda,  1988]  that  the  high  latitude  ionosphere  and 
magnetosphere  can  be  characterized  as  a  highly  turbulent  and  structured 
plasma  containing  density,  electric  field,  and  magnetic  field  fluctuations 
with  scale  sizes  ranging  from  approximately  hundreds  of  kilometers  to 
centimeters. 

Several  different  source  mechanisms  [Kintner  and  Seyler,  1985;  Temerin 
and  Kintner,  1988;  Tsunoda,  1988]  have  been  proposed  to  account  for  high 
latitude  magnetospheric  and  ionospheric  plasma  turbulence,  e.g. ,  electric 
fields,  particle  precipitation,  and  plasma  instabilities.  The  relative 
importance  of  each  mechanism  depends  on  several  factors,  e.g.,  geomagnetic 
latitude,  seasonal  effects,  and  scale  size  perpendicular  to  the  geomagnetic 
field.  Keskinen  and  Ossakow  [1982,  1983]  have  proposed  the  E  x  B  drift 
instability  to  be  a  source  of  high  latitude  ionospheric  plasma  density  and 
electric  field  fluctuations  for  scale  sizes  from  tens  of  kilometers  to 
meters.  This  interchange-like  instability,  a  plasma  analogue  of  the  well- 
known  Rayleigh-Taylor  instability  which  results  when  a  heavy  fluid  is 
supported  against  gravity  by  a  lighter  fluid,  is  driven  by  both  a  plasma 
density  gradient  and  an  electric  field  perpendicular  to  the  geomagnetic 
field.  The  electric  field  will  result  in  a  relative  velocity  between 
convecting  high  latitude  ionospheric  plasma  ions  and  the  neutral 
thermosphere.  The  E  x  B  drift  instability  has  been  applied  to  the  small 
scale  dynamics,  stability  and  evolution  of  convecting  large  scale 
ionospheric  density  structures  recently  observed  in  the  high  latitude 
ionosphere  [Vickrey  et  al.,  1980;  Muldrew  and  Vickrey  1982;  Buchau  et  al., 
1983,  1985;  Weber  et  al.,  1984,  1986;  Foster  and  Doupnik,  1984;  de  la 
Beaujardiere  et  al.,  1985] .  Several  studies  [Sojka  and  Schunk,  1986; 
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Schunk  and  Sojka,  1987;  Anderson  et  al.,  1987]  have  modeled  the  large 
scale,  global  features  of  these  plasma  density  enhancements  and  structures. 
Vickrey  and  Kelley  [1982]  have  considered  the  effects  of  a  conducting  E 
layer  on  the  evolution  of  convecting  large  scale  ionospheric  density 
structures.  At  smaller  scales,  Huba  et  al  [1983]  have  studied  the  linear 
theory  of  the  E  x  B  instability  with  an  inhomogeneous  electric  field. 
Mitchell  et  al.  [1985]  have  investigated  the  nonlinear  evolution  of  the  E  x 
B  instability  in  the  high  latitude  ionosphere  with  magnetospheric  coupling. 
They  find  that  the  primary  effect  of  the  magnetosphere  is  to  incorporate 
inertia  into  the  development  of  the  E  x  B  instability.  In  addition  they 
find  that  plasma  interchange-like  instabilities  in  the  high  latitude 
ionosphere  develop  in  a  fundamentally  different  manner  when  magnetospheric 
coupling  is  included  as  opposed  to  the  case  when  magnetospheric  coupling  is 
absent.  Chaturvedi  and  Huba  [1987]  have  added  three-dimensional  effects  to 
the  linear  theory  of  the  E  x  B  instability  in  the  high  latitude  ionosphere. 
It  is  now  generally  believed  [Temerin  and  Kintner,  1988;  Tsunoda,  1988] 
that  plasma  interchange-like  instabilities  are  a  major  source  of  density 
and  electric  field  fluctuations  in  the  high  latitude  ionosphere  in  the 
scale  size  regime  of  a  few  tens  of  kilometers  to  tens  of  meters 
perpendicular  to  the  geomagnetic  field. 

However,  the  spectrum  of  density  or  electric  field  fluctuations 
associated  with  the  nonlinear  evolution  of  the  E  x  B  instability  in  the 
coupled  magnetosphere-ionosphere  system  has  not  been  analytically 
investigated  in  detail.  In  this  study  we  derive  constraints  on  and  scaling 
laws  associated  with  the  nonlinear  spectrum  of  the  E  x  B  interchange 
instability  as  it  may  occur  in  the  high  latitude  ionosphere-magnetosphere 
coupled  system.  These  constraints  are  derived  by  analysis  of  the 
fundamental  nonlinear  equations  governing  the  E  x  B  instability  in  Sec.  2. 
We  discuss  and  summarize  our  results  in  Sec.  3. 

2.  MODEL  EQUATIONS  AND  ANALYSIS 

We  consider  the  high  latitude  ionospheric  plasma  to  be  a  weakly  ionized 
low  3  =  8ilng(Tg  +  Tj)/Bq^  plasma  where  n^  is  the  electron  plasma  density, 
Te(Ti)  the  electron  (ion)  temperature,  and  Bp  the  geomagnetic  field.  For 
our  coordinate  system  we  take  the  ambient  density  gradient  to  be  in  the  x- 
direction,  i.e.,  n  =  n^Cx),  the  ambient  electric  field  to  be  in  the  y- 
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direction,  i.e.,  E  =  EqZ*  geomagnetic  field  to  be  in  the  z- 
direction  B  =  BqZ.  The  equations  used  in  this  analysis  are  continuity, 
momentum,  and  charge  conservation: 


an 

^  +  V-n  V  =  0 
dt  ot-a 


(1) 


0  =  -  — Te  +  c  X  b1  -  (T  /m  )(7n/n)  - 
m  -e  -)  '  e  e 


V  .(V  -  V.) 

ei'-e  -1 


(2) 


(Ir  »  -  |r(l  *  <^'hi  »  b)  -  <T,/n,.)(Vn/n)  -  -  V^)  -  (3) 


V-J  = 


=  0 


(A) 


Here,  n^j^  is  the  electron  (a=e)  or  ion  (a=i)  density,  is  the  electron  or 
ion  velocity,  E  is  the  electric  field,  e  is  the  electron  charge,  m^  is  the 
electron  or  ion  mass,  c  is  the  speed  of  light,  is  the  electron  or  ion 
temperature,  '’ei^^ie^  electron-ion  (ion-electron)  collision 
frequency,  and  is  the  ion-neutral  collision  frequency. 

We  now  perturb  Eq.  (l)-(4)  about  an  equilibrium  and  let  n^^  =  uq  +  Sn,  E 
=  Eo  -  WS()>,  and  +  6V  with  Sn,  5E,  5V  «  exp  [i(kjjX  +  kyy)  -  (ot ] 
where  (k^  +  ky)^'^^  L  >  1  and  =  nQ^(anQ/ax).  After  transforming  to  a 
reference  frame  moving  with  velocity  =  (cEp  x  B)/B^  =  (cEq/B)x,  Eqs.  (2) 
and  (3)  yield. 


V  =  V  + 
-e  -eo 


5V 


-e 
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c  1_  f  3_ 
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where  &?  =  6«t>  -  (Tg/e)  (  Sh/Hq)  ,  6<J;  =  5<|>  +  (T^/e)(  Sn/ng) ,  9g(9^)  is  the 
electron  (ion)  gyrof requency ,  and  Vg(vj)  is  the  electron  (ion)  thermal 
velocity.  Mitchell  et  al  [1985]  and  Keskinen  et  al  [1988]  have  shown  that 
electric  fields  generated  by  the  E  x  B  instability  in  the  high  latitude 
ionosphere  can  map  far  into  the  magnetosphere  on  time  scales  comparable  to 
the  E  X  B  instability  growth  time.  As  a  result,  a  magnetic-field-line 
integrated  model  of  the  E  x  B  instability  is  appropriate  for  application  to 
the  high  latitude,  near-earth  space  plasma  [Mitchell  et  al.,  1985].  Using 
Eqs.  (1)  and  (4)  one  can  then  write  [Mitchell  et  al.,  1985], 
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V  =  -  |7$  X  2 


(9) 
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a 
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with  =  B/(4iinmj)^72  Alfvfen  speed,  D  =  ( Vg )  [  1+Tg/Tj  ] ,  and  E  = 
with  $  =  ((i  -  (Tg/e)  In  n.  Within  the  context  of  this  model,  Eq.  (8) 
then  defines  an  effective  frequency  v  =  Ep/C^.  For  fluctuations  with 
characteristic  frequency  a/at  ~  y,  then  v/y  is  essentially  the  flux-tube- 
integrated  ratio  of  the  Pedersen  current  in  the  ionosphere  and  polarization 
current  in  the  magnetosphere.  The  quantity  v/ y  <  1  implying  magnetospheric 
control  for  small  Pedersen  conductivities,  large  inertial  capacitances,  and 
large  fluctuation  growth  rates.  On  the  other  hand,  v/y  >  1  implying 
ionospheric  control  for  large  Pedersen  conductances,  small  inertial 
capacitances,  and  small  fluctuation  growth  rates.  In  order  to  solve  Eq. 
(7)-(8)  one  must  specify  Op  and  along  a  specified  flux  tube.  We  make 

the  assumption  that  the  magnetosphere,  to  lowest  order,  is  uniform  and 
incompressible.  As  a  result,  a  continuity  equation  for  the  plasma  density 
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in  the  magnetosphere  is  not  necessary.  We  further  take  the  Pedersen 
current  Ep  'oN  and  the  inertial  capacitance  C^,  “  N  with  N  =  J  n  dz  the 
integrated  plasma  density.  We  make  the  assumption  that  the  magnetosphere 
is  uniform  and  incompressible.  Eqs.  (7)-(9)  can  be  further  simplified 
giving: 

§7  +  7-NV  =  DW^N  (10) 


W-  +  n(|^  +  V-W)W$  +  (Ti/e)vWNj  =  0  (11) 


V  =  -  I  7$  X  z  (12) 

Eq.  (10-(12)  are  now  a  closed  system  in  N  and 

Eq.  (lO)-(ll)  have  been  used,  in  the  limit  T^,  Tg  =  0,  to  study  the 
nonlinear  evolution  of  the  interchange  E  x  B  instability  in  the  high 
latitude  ionosphere  by  Mitchell  et  al.  11985).  They  find  that  the  E  x  B 
in'itabi  1 3  develops  in  a  fundamentally  different  manner  depending  on 
whether  the  collisional  (v/y  >  1)  or  inertial  (v/y  <  1)  limit  is  taken. 

Here  y  is  the  linear  growth  rate  of  the  E  x  B  instability.  We  now  proceed 
to  derive  scaling  laws  associated  with  spectrum  of  the  interchange  E  x  B 
instability  in  the  inertial  limit.  Expanding  N  =  n^fx)  +  5n(x,y)  and  ?  = 
<))q  +  6<)>(x,  y)  we  find 

+  5-  z  X  y5(t>-Vn  +  I  z  X  W6())- V6n  =  DV^6n  (13) 

O  t  D  O  D 

Multiplying  Eq.  (13)  by  5n  and  integrating  over  all  x  and  y  with  d^x  =  dxdy 
we  find 
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j  ^  &n  z  X  VS(|c  7n^ 


^  [d^x  5n  z  X  75<<i*  75n  -  D  fd^x  5n7^5n  =  0  (1^) 


In  Eq.  (14),  the  second  term  on  the  left-hand  side,  resembling  a  source 
term,  is  proportional  to  the  initial  density  gradient  while  the  fourth 
term,  resembling  a  sink,  is  proportional  to  the  cross-field  diffusion.  The 
third  term,  which  is  nonlinear  and  represents  a  mode-mode  coupling  among 
the  fluctuations,  vanishes  upon  averaging  over  d^x  since 


J' 


6n  z  X  76<|>'  Sn 


■ 

d  X  z 


X  7S<f7(&n)‘ 


=  j  Jd^x  7-  (Sn)^z  X  VS^]  =  0 


since  (6n)2  -♦  0  as  x,  y  -♦  ”  for  a  finite  sized  ionosphere  where  7-((6n)2  z 

X  75(frJ  =  z  X  75<fr-7(5n)^.  We  Fourier  expand  6n  and  as  follows: 


'Sn(x,y)' 

d  k 

-6n, 

.54>(x,y). 

« 

.&<l>k' 

exp  ik'X 


with  k  =  kjjX  +  ky£.  As  a  result,  Eq.  (14)  can  be  written 
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.d|, 


k^i, 


=  0 


(15) 


with  I, 
k 


Eq.  (11)  gives 


S^i. 


=  n  k^( iw  -  v)  ^ [ 
.0  J  L 


V  (ik-Eo  + 


(Ti/e)k2)] 


6n, 


(16) 


which  can  be  inserted  into  Eq.  (15).  We  not  that  Eq.  (15)  using  Eq.  (16) 
can  be  analyzed  either  in  the  collisional  (v  >  w)  limit,  as  was  done  by 
Keskinen  and  Ossakow  [1981],  or  the  inertial  (v  <  u)  limit.  In  the 
inertial  limit  we  have 

?  It  ■  Fi^Vk 

with 

\  =  (h  *  lo  -  ’  ^"o)]  -  Dk^  =  (  vcEpCOS  e/BL)^'^^-  Dk^ 

with  L”1  =  ( I/Oq) ( SOq/Sx)  and  0  is  the  angle  defined  by  k  and  Eq. 

Eq.  (17)  can  be  written  in  polar  coordinates  (k^  =  kj^^  +  ky2,  e  =  tan“l 
kj{/ky)  assuming  steady  state  conditions  in  the  following  manner: 


C 

max 

dk  k 

d0  Y.  „  l(k,0)  = 

dk  k 

•'k  . 

0  ^ 

k 

min 


r2n 

0 


de  Y,^^^i(k,0)  (18) 


where  =  ( vcos0cEq/BL)1^2  jg  growth  rate  and  =  -Dk^  is  the 
damping  rate.  As  a  result,  the  spectral  power  generated  in  the  unstable 
range  of  wave  numbers  between  kn,iQ  and  k^.,  the  critical  wave  number,  is 
balanced  by  the  power  dissipated  in  the  stable  range  of  wave  numbers 
between  k^,  and  k^gjj  under  steady  state  conditions. 
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In  the  following  we  assume  a  general  form  for  the  density  spectrum 
l(k.,0)  and  use  Eq.  (18)  to  find  constraints  on  the  parameters  used  to 
specify  l(k,0).  To  our  knowledge,  the  steady  state  spectrum  l(k,0)  of  the 
E  X  B  instability  in  the  high  latitude  ionosphere,  in  the  inertial  regime, 
has  not  been  studied  in  detail.  Keskinen  and  Ossakow  (1982,  19831  have 
studied  aspects  of  the  spectrum  in  the  collisional  regime.  However,  Hassam 
et  al.  [1986]  have  investigated  the  nonlinear  spectrum  of  the  gravitational 
Rayleigh-Taylor  instability,  a  closely  related  instability  to  the  E  x  B 
instability,  in  both  the  inertial  and  collisional  regimes.  Hassam  ec  al. 
[1986]  show  that  the  spectrum  in  the  inertial  regime  is  more  isotropic  than 
in  the  collisional  regime.  take,  consistent  with  the  computations  of 

Hassam  et  al.  [1986],  the  spectrum  to  be  of  the  form 

l(k,0)  =  lof(0)  (1  +  (k/ko)2)-(n+l)/2 

with  f(0)  =  cos''’0  +  h  sin"*©,  q  <  1,  and  m  >  0.  Substituting  l(k,0)  into 

Eq.  (18)  we  find 


3-n 

I'm  +  2'] 

1  2  J 

r 

fm  +  2^ 

1  2  J 

1  1/2 
(vcE  /BL) 

n-1  j.| 

fm  +  1) 

1  2  J 

r 

m  +  I) 

1  2  J 

1  Dk^ 

'  max 

where  r(q)  is  the  gamma  function  defined  by 


r(q)  =  J  dy  y'^  ^e  ^ 
0 


(19) 


and  we  have  assumed  n  >  1.  From  Eq.  (19)  we  see  that  1  <  n  <  3.  We  note 
that  the  scaling  in  the  inertial  limit  for  k^  as  given  by  Eq.  (19)  is 
different  from  the  collisional  limit  (Keskinen  and  Ossakow,  1981]  where 
(kg/k^gj^)®  (cEq/BLD  .  In  order  to  verify  the  spectral  scaling 
laws  outlined  in  this  paper,  density  power  spectra  in  or  near  converting 
patches  and  plasma  enhancements  in  the  high  latitude  ionosphere  must  be 
measured.  The  spectral  index,  outer  scale  k^,  and  inner  scale  kj^gj^ 
wavenumbers  would  then  need  to  be  computed.  The  scaling  of  l<o/*^max 
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the  ambient  electric  field  Eq  in  the  neutral  frame,  say,  could  then  be 
investigated.  Taking  n  =  2,  m  =  2,  (vBL/cEq)  =  0.2  [Mitchell  et  al., 
1985],  cEq/B  =  500  m/s,  D  =  1  m^/S,  L  =  10km  and  2ii/k^gjj  =  10  m  [R.T. 
Tsunoda,  1989]  we  find  from  Eq.  (19)  2n/kQ  =  600  m  which  is  consistent  with 
recent  observations  (R.T.  Tsunoda,  1989]. 

The  gravitational  Rayleigh-Taylor  instability,  thought  to  be 
responsible  for  rising  bubble  and  plumelike  phenomena  in  the  equatorial 
ionosphere,  [Ossakow,  1981;  Kelley  and  McClure,  1981]  is  analogous  to  the 
interchange  Ex  B  instability.  The  Rayleigh-Taylor  instability  can  be 
described  by  equations  identical  to  Eq.  (l)-(4)  except  that  the  driving 
electric  field  Eq  =  Eq^  is  replaced  by  gravity  g  =  -gx.  The  formulae  for 
Rq  in  Eq.  (19)  which  are  applicable  to  interchange  instability  due  to  the 
gravitational  Rayleigh-Taylor  instability  may  then  be  computed  simply  by 
taking  Eq  -»  Bg/cvj^Q.  Recently,  experimental  observations  gathered  during 
the  PLUMEX  rocket  campaign  (Rino  et  al,  1981;  Kelley  et  al.,  1982]  have 
been  further  analyzed  (R.C.  Livingston,  private  communication].  These  data 
indicate  that  the  power  spectrum  of  density  fluctuations  in  equatorial 
Spread-F  usually  contains  a  break  at  a  wavelength  Xq  which  is  altitude 
dependent.  The  power  spectrum  is  steeper  (shallower)  for  wavelengths 
greater  (less)  than  the  break  wavelength.  Typical  power  spectra  Pj^  a  k~'^ 
have  spectral  indices  n  =  1.5  -  3  and  are  computed  in  the  altitude  interval 
of  approximately  275-500  km.  Furthermore,  the  data  indicate  [R.C. 
Livingston  private  communication]  that  Xg  is  directly  proportional  to 
altitude,  i.e.,  the  break  wavelength  Xg  is  large  (small)  at  high  (low) 
altitudes  with  Xg  =  0.2  -  1  km  in  the  altitude  interval  z  =  275-500  km.  In 
addition,  it  was  found  [R.  C.  Livingston,  private  communication]  that  the 
altitude  dependence  given  by  the  ion-neutral  collision  frequency, 

provides  a  good  fit  to  the  data,  i.e.,  Xg  =  or  kg  with  kg  = 

2n/Xg.  Using  Eq.  (19)  and  the  previous  results  of  Keskinen  and  Ossakow 
[1981]  we  can  then  write  for  the  gravitational  Rayleigh-Taylor  instability, 
('^o^*^max)"~^  “  inertial  limit,  i.e.,  <  (g/L)^^^ 

while  (ItQ/k^ax)'’”^  “  S/'’in^^*^max  collisional  limit  >  (g/L)^^^. 

These  scalings  can  be  rewritten  as  (l^o^^^max^"^  ^  “  IC/ ''"D^ *^max ^  where  yg  is 
the  growth  rate  and  YD(k^gx)  is  the  damping  rate  for  mode  k  evaluated  at  k 
=  k^a^.  As  a  result  kj/k^  where  kj(k^)  is  the  spectral  break 

wavenumber  in  the  inertial  (collisional)  limit  and  yJ  =  (g/L)^^^  is  the 
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growth  rate  in  the  inertial  limit  with  growth  rate  in  the 

collisional  limit.  Since  yJ  <  y^»  typically,  kp/k,^  <  1,  i.e.,  the  spectral 

break  wavelength  XQ«2n/kQ  is  larger  in  the  inertial  regime  (high  altitudes) 

than  in  the  collisional  regime  (low  altitudes)  in  agreement  with  the 
observations  of  Livingston  et  al.  Furthermore,  the  scaling  of  the  spectral 
break  wavenumber  with  altitude  kj/kg  yJ/y^  “  is  also  consistent  with 
the  PLUMEX  data  [R.C.  Livingston,  private  communication]. 

3.  SUMMARY 

We  have  derived  constraints  and  scaling  relationships  for  the  power 
spectrum  of  the  inertial  E  x  B  drift  instability  in  the  weakly  nonlinear 
regime  in  the  high  latitude  ionosphere  by  analyzing  conservation  laws 
implied  by  the  fundamental  plasma  fluid  equations.  By  assuming  l(k,0)  ® 
f(0)(l  +  (k/kQ)^)"^'’'^^^'^^  where  I  (k,0)  is  the  two-dimensional  density 

power  spectrum  of  the  E  x  B  instability,  cos  0  =  k'Eg,  Eq  is  the  ambient 

electric  field,  f(0)  =  cos'"0  +  q  sin‘''0,  m,  q  constants,  and  k^  is  the  outer 
scale  we  find  that  1  <  n  <  3.  Furthermore,  we  derive  the  scaling  of  the 
ratio  of  outer  (kg)  and  inner  (k^jgj^)  scales  kQ/k,„gj(.  «  ( vcEQ/BLD^k^3x)^'^^^'^~ 
with  k^gjj  the  inner  scale.  Using  typical  values  of  v,  Eq,  L,  D,  and 
kmax  find  values  for  kg  and  n  which  are  consistent  with  observations. 
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ATTN 

OPR/ALVA  T.  STAIR 

OICY 

ATTN 

PHD/ JURGEN  BUCHAU 

OICY 

ATTN 

PHD/ JOHN  P.  MULLEN 

AF  WEAPONS  LABORATORY 
KIRTLAND  AFB,  NM  87117 
OICY  ATTN  SUL 

OICY  ATTN  CA/ARTHUR  H.  GUENTHER 

AFTAC 

PATRICK  AFB,  FL  32925 
OICY  ATTN  TN 

WRIGHT  AERONAUTICAL  LABORATORIES 
WRIGHT-PATTERSON  AFB,  OH  45433-6543 
OICY  ATTN  AAAI /WADE  HUNT 

OICY  ATTN  AAAI /ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  AND  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  DC  20330 
OICY  ATTN  AFRDQ 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731-5000 
OICY  ATTN  J.  DEAS 

ESD/SCD-4 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
OICY  ATTN  NI CD/ LIBRARY 

OICY  ATTN  ETDP/B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFIN  AFB,  NY  13441 

OICY  ATTN  DOC  LIBRARY/TSLD 

OICY  ATTN  OCSE/V.  COYNE 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 
OICY  ATTN  XPFS 

SAMSO/MN 

NORTON  AFB,  CA  02409 
(MINUTEMAN) 

OICY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EEP/A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY,  ROOM  G-042 
WASHINGTON,  DC  20545 

OICY  ATTN  DOC  CON  FOR 

A.  LABOWITZ 

DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR 

D.  SHERWOOD 


EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR 
J.  BREEDLOVE 
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UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 


OICY 

ATTN 

DOC  CON  FOR 
TECH  INFO  DEPT 

OICY 

ATTN 

DOC  CON  FOR 
L-389/R.  OTT 

OICY 

ATTN 

DOC  CON  FOR 
L-31/R.  HAGER 

LOS  ALAMOS  NATIONAL  LABORATORY 
P.O.  BOX  1663 

LOS  ALAMOS,  NM 

87545 

OICY 

ATTN 

J.  WOLCOTT 

OICY 

ATTN 

E.  JONES 

OICY 

ATTN 

J.  MALIK 

OICY 

ATTN 

R.  JEFFRIES 

OICY 

ATTN 

J.  ZINN 

OICY 

ATTN 

D.  WESTERVELT 

OICY 

ATTN 

D.  SAPPENFIELD 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 
AEROSPACE  CORPORATION 


1.  BOX  92957 
;  ANGELES,  CA 

90009 

OICY 

ATTN 

I.  GARFUNKEL 

OICY 

ATTN 

T.  SALMI 

OICY 

ATTN 

S.  BOWER 

OICY 

ATTN 

D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOCIATION,  INC. 
1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

OICY  ATTN  L.  SLOAN 

OICY  ATTN  R.  THOMPSON 


LOS  ALAMOS  NATIONAL  LABORATORY 
MS  D438 

LOS  ALAMOS,  NM  87545 


OICY 

ATTN 

S.P.  GARY 

OICY 

ATTN 

J.  BOROVSKY 

SANDIA  LABORATORIES 

P.O.  BOX  5800 

ALBUQUERQUE,  NM 

87115 

OICY 

ATTN 

W.  BROWN 

OICY 

ATTN 

A.  THORNBROUGH 

OICY 

ATTN 

T.  WRIGHT 

OICY 

ATTN 

D.  DAHLGREN 

OICY 

ATTN 

3141 

OICY 

ATTN 

SPACE  PROJ  DIV 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 


OICY 

ATTN 

J.  WORKMAN 

OICY 

ATTN 

C.  PRETTIE 

OICY 

ATTN 

S.  BRECHT 

BOEING  COMPANY, 
P.O.  BOX  3707 

THE 

SEATTLE, 

WA  98124 

OICY 

ATTN 

G.  KEISTER 

OICY 

ATTN 

D.  MURRAY 

OICY 

ATTN 

G.  HALL 

OICY 

ATTN 

J.  KENNEY 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  B.  MURPHEY 

OICY  ATTN  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  DC  20545 

OICY  ATTN  DR.  YO  SONG 


CHARLES  STARK  DRAPER  LABORATORY 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  92139 

OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
22300  COMSAT  DRIVE 
CLARKSBURG,  MD  20871 
OICY  ATTN  G.  HYDE 


NATL.  OCEANIC  &  ATMOSPHERIC 
ADMINISTRATION 
ENVIRONMENTAL  RESEARCH  LABS 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

OICY  ATTN  R.  GRUBB 
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CORNELL  UNIVERSITY 

DEPT  OF  ELECTRICAL  ENGINEERING 

ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

OICY  ATTN  H.  LOGSTON 

OICY  ATTN  SECURITY/ 

(PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  WILSHIRE  BLVD. 

SANTA  MONICA,  CA  90401 

OICY  ATTN  C.G.  GABBARD 

OICY  ATTN  R.  LELEVIER 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

OICY  ATTN  SECURITY  OFFICER 

OICY  ATTN  T.N.  DAVIS 

OICY  ATTN  NEAL  BROWN 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP- 
EASTERN  DIVISION 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

OICY  ATTN  DAN  MCCLELLAND 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 

OICY  ATTN  ERNEST  BAUER 

OICY  ATTN  HANS  WOLFARD 

INTL  TELL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 


JAYCOR 

P.O.  BOX  85154 
11011  TORRE YANA  ROAD 
SAN  DIEGO,  CA  92138 

OICY  ATTN  N.T.  GLADD 

OICY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOC  LIBRARIAN 

OICY  ATTN  THOMAS  POTEMRA 

OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORPORATION 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET 
(P.O.  DRAWER  00) 

SANTA  BARBARA,  CA  93102 


OICY 

ATTN 

DASIAC 

OICY 

ATTN 

WARREN  S.  KNAPP 

OICY 

ATTN 

WILLIAM  MCNAMARA 

OICY 

ATTN 

B.  GAMBILL 

LINKABIT  CORPORATION 

10453  ROSELLE 

SAN  DIEGO,  CA  92121 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT/ 

DEPT  52-12 

OICY  ATTN  W.L.  IMHOF/ 

DEPT.  52-12 

OICY  ATTN  J.B.  CLADIS/ 

DEPT.  52-12 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  02647 
OICY  ATTN  W.  OLSON 

OICY  ATTN  R.W.  HALPRIN 

OICY  ATTN  TECHNICAL 

LIBRARY  SERVICES 
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MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 


SANTA  BARBARA, 

CA  03101 

OICY 

ATTN 

P.  FISCHER 

OICY 

ATTN 

W.F.  CREVIER 

OICY 

ATTN 

STEVEN  L.  GUTSCHE 

OICY 

ATTN 

R.  BOGUSCH 

OICY 

ATTN 

R.  HENDRICK 

OICY 

ATTN 

RALPH  KILB 

OICY 

ATTN 

DAVE  SOWLE 

OICY 

ATTN 

F.  FA JEN 

OICY 

ATTN 

M.  SCHEIBE 

OICY 

ATTN 

CONRAD  L.  LONGMIRE 

OICY 

ATTN 

B.  WHITE 

OICY 

ATTN 

R.  STAGAT 

OICY 

ATTN 

D.  KNEPP 

OICY 

ATTN 

C.  RINO 

MISSION 

RESEARCH  CORPORATION 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NM 

87106 

OICY 

ATTN 

R.  STELLINGWERF 

OICY 

ATTN 

M.  ALME 

OICY 

ATTN 

L.  WRIGHT 

MITRE  CORPORATION 
WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

OICY  ATTN  W.  HALL 

OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
123A0  SANTA  MONICA  BLVD 
LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
UNIVERSITY  PARK,  PA  16802 
01  CY  ATTN  IONOSPHERIC 

RESEARCH  LAB 

PHOTOMETRICS,  INC. 

A  ARROW  DRIVE 
WOBURN,  MA  01801 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  04610 


R  &  D  ASSOCIATES 


P.O.  BOX  9695 
MARINA  DEL  REY, 

CA  90291 

OICY 

ATTN 

FORREST  GILMORE 

OICY 

ATTN 

W.B.  WRIGHT,  JR 

OICY 

ATTN 

W.J.  KARZAS 

OICY 

ATTN 

H.  ORY 

OICY 

ATTN 

C.  MACDONALD 

OICY 

ATTN 

BRIAN  LAMB 

OICY 

ATTN 

MORGAN  GROVER 

RAYTHEON  CORPORATION 
528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42ND  STREET 
NEW  YORK,  NY  10036 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
10260  CAMPUS  POINT  DRIVE 
SAN  DIEGO,  CA  92121-1522 
OICY  ATTN  L.M.  LINSON 

OICY  ATTN  D.A,  HAMLIN 

OICY  ATTN  E.  FRIEMAN 

OICY  ATTN  E.A.  STRAKER 

OICY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
1710  GOODRIDCE  DRIVE 
MCLEAN,  VA  22102 

OICY  ATTN  J.  COCKAYNE 

OICY  ATTN  E.  HYMAN 
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SRI  INTERNATIONAL 
333  RAVENSUOOD  AVENUE 
MENLO  PARK,  CA  94025 


OICY 

ATTN 

J.  CASPER 

OICY 

ATTN 

DONALD  NEILSON 

OICY 

ATTN 

ALAN  BURNS 

OICY 

ATTN 

G.  SMITH 

OICY 

ATTN 

R.  TSUNODA 

OICY 

ATTN 

D.A.  JOHNSON 

OICY 

ATTN 

U.G.  CHESNUT 

OICY 

ATTN 

C.L.  RINO 

OICY 

ATTN 

WALTER  JAYE 

OICY 

ATTN 

J.  VICKREY 

OICY 

ATTN 

R.L.  LEADABRAND 

OICY 

ATTN 

G.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

R.  LIVINGSTON 

OICY 

ATTN 

V.  GONZALES 

OICY 

AT  I'N 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

OICY  ATTN  R.K.  PLEBUCH 

OICY  ATTN  D.  DEE 

OICY  ATTN  D.  STOCKWELL/ 

SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 


BURLINGTON 

1,  MA 

01803 

OICY 

ATTN 

W. 

REIDY 

OICY 

ATTN 

J. 

CARPENTER 

OICY 

ATTN 

C. 

HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

OICY  ATTN  N.  ZABUSKY 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE  (UNLESS 
OTHERWISE  NOTED) 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  DC  20375-5000 

DR.  S.  OSSAKOW,  CODE  4700  (2  CYS) 
CODE  4701 

DR.  J.  HUBA,  CODE  4780  (2  CYS) 

DR.  H.  GURSKY,  CODE  4100 
DR.  J.M.  GOODMAN,  CODE  4180 
DR.  P.  RODRIGUEZ,  CODE  4750 
DR.  P.  MANGE,  CODE  1004 
DR.  R.  MEIER,  CODE  4140 
CODE  2628  (22  CYS) 

CODE  1220 

A.F.  GEOPHYSICAL  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01731 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  W.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 
DR.  J.F.  RICH 
DR.  N.  MAYNARD 
DR.  D.N.  ANDERSON 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  02215 
DR.  J.  AARONS 
DR.  M.  MENDILLO 
DR.  J.M.  FORBES 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
1801  N.  BEAUREGARD  STREET 
ARLINGTON,  VA  22311 
DR.  E.  BAUER 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
CAMBRIDGE,  MA  02139 
LIBRARY,  NW16-262 
DR.  T.  CHANG 
DR.  R.  LINDZEN 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 

DR.  R.F.  PFAFF,  CODE  696 

DR.  R.F.  BENSON 

DR.  K.  MAEDA 

DR.  S.  CURTIS 

DR.  M.  DUBIN 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  95152 

MR.  R.  ROSE,  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND 
ENVIRONMENTAL  LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  K.  DAVIES 
DR.  R.F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 
DR.  C.  ROBERSON 

LABORATORY  FOR  PLASMA  AND 
FUSION  ENERGIES  STUDIES 
UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20742 
JEAN  VAFYAN  HELLMAN, 

REFERENCE  LIBRARIAN 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  E.  KLEVANS 
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PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 

SAIC 

10260  CAMPUS  POINT  DRIVE 
SAN  DIEGO,  CA  92121-1522 
DR.  L.  LINSON/MS  33 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  04025 
DR.  R.  TSUNODA 
DR.  WALTER  CHESNUT 
DR.  J.  VICKREY 
DR.  R.  LIVINGSTON 

STANFORD  UNIVERSITY 
STANFORD,  CA  04305 
DR.  P.M.  BANKS 
DR.  R.  HELLIWELL 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AL  99701 
DR.  L.C.  LEE 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UT  84322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 
DR.  B.  FEJER 

UNIVERSITY  OF  CALIFORNIA 
F.OS  ALAMOS  NATIONAL  LABORATORY 
EES  DIVISION 
LOS  ALAMOS,  NM  87545 

DR.  M.  PONGRATZ,  EES-DOT 
DR.  D.  SIMONS,  ESS-7,  MS-D466 
DR.  S.P.  GARY,  ESS-8 
DENNIS  RIGGIN,  ATMOS  SCI  GROUP 


UNIVERSITY  OF  ILLINOIS 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
1406  V.  GREEN  STREET 
URBANA,  IL  61801 
DR.  ERHAN  KUDEKI 

UNIVERSITY  OF  CALIFORNIA, 

LOS  ANGELES 
405  HILLGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR.  C.  KENNEL 


UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 
DR.  R.  GREENWALD 
DR.  C.  MENG 
DR.  T.  POTEMRA 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 

UNIVERSITY  OF  TEXAS  AT  DALLAS 
CENTER  FOR  SPACE  SCIENCES 
P.O.  BOX  688 
RICHARDSON,  TX  75080 
DR.  R.  HEELIS 
DR.  W.  HANSON 
DR.  J.P.  McCLURE 

DIRECTOR  OF  RESEARCH 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 
(2  CYS) 

UNIVERSITY  OF  SAGAR 
DEPARTMENT  OF  PHYSICS 
SAGAR-4 70003 
(M.P.)  INDIA 

DR.  M.S.  TIWARI 


26 


